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Abstract 

Pair creation of scalar particles in a boost-invariant electric field which is confined in the 
forward light cone is studied. We present the proper-time evolution of momentum distribu- 
tions of created particles, which preserve the boost invariance of the background field. The 
two-particle correlation of the created particles is also calculated. We find that long-range 
rapidity correlations may arise from the Schwinger mechanism in the boost-invariant electric 
field. 

1 Introduction 

Boost invariance in a longitudinal beam direction is a key ingredient in theoretical descriptions 
of ultrarelativistic heavy-ion collisions. In 1983, Bjorken pQ introduced a boost-invariant initial 
condition for Landau's hydrodynamic model [2j. Because the hydrodynamic equations as well 
respect the boost symmetry, the following space-time evolution of the system is also boost- 
invariant, which leads to the central plateau in the particle rapidity distributions. Owing to 
this symmetry, the fluid evolution is described by the fluid proper time r = \/t 2 — z 2 and the 
transverse coordinates xj_ and does not depend on the space-time rapidity r\ = \ In |i| , where 
the longitudinal beam direction is taken along the z axis. 

Before the hydrodynamic evolution, the system of ultrarelativistic heavy-ion collisions is also 
assumed to hold the boost symmetry, because shortly after the collision, the nuclei are highly 
Lorentz-contr acted pancakes receding from the collision point at nearly the speed of light, so 
that there is no scale in the longitudinal direction. For example the flux-tube model (the Low- 
Nussinov model O d]), which assumes the generation of longitudinal color-electric fields just 
after high-energy particle collisions, comprehends the boost invariance. Also in the framework 
of the color glass condensate, the formation of boost-invariant longitudinal electric fields and 
the longitudinal magnetic fields as well is predicted, of which the state is called glasma f5J. Such 
electric fields would depend on r and xj_ but do not depend on the space-time rapidity 77. 

Starting from the boost-invariant classical fields, how does the system evolve into a boost- 
invariant and locally thermalized quark-gluon plasma? One possible approach to this problem 
is solving the classical Yang-Mills equations and investigating the dynamics of the classical 
electromagnetic fields Because of the expansion of the system, the electromagnetic 

fields decay in time. To understand the formation process of the quark-gluon plasma, however, 
the notion of particle or quantized fields is necessary. In particular, the classical Yang-Mills 
equations cannot describe the quark production. If one introduces charged quantum fields 
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interacting with the classical electric field into the system evolution, the system becomes unstable 
against particle production via the Schwinger mechanism [8]. 

The Schwinger mechanism, nonperturbative particle pair creation from a classical field, has 
been studied as a mechanism of matter formation in heavy-ion collisions [9l \10\ lllj . However, 
the original Schwinger mechanism treats spatially uniform electric fields. In contrast, electric 
fields expected to be generated in the early stage of heavy-ion collisions exist only between two 
nuclei receding from each other at a velocity close to the speed of light. In an ideal situation 
where the two nuclei run at exactly the speed of light, the electric fields span only inside the 
forward light cone. 

In the case of spatially homogeneous fields, the semiclassical tunneling description of pair cre- 
ation implies that particles and antiparticles are created with vanishing longitudinal momentum. 
This has been assured by field-theoretical calculations [12]. Clearly, these particle-antiparticle 
pairs break the boost invariance of the background field by introducing one specific frame, 
namely their center of mass frame. How is this situation modified if an electric field exists only 
inside the light cone? 

In this paper, we investigate the Schwinger mechanism under an electric field which localizes 
in the forward light cone, focusing on its real-time dynamics and the Lorentz-boost invariance of 
the system. To get the real-time description, we introduce time-dependent particle definitions 
\12\ 113] . The canonical quantization is done on an equal-r line instead of an equal-time line in 
order to give a description preserving the boost invariance. The proper quantization procedure in 
the boost-invariant curvilinear coordinates [13] provides a boost-invariant description of the pair 
creation. To concentrate on revealing the boost-invariant dynamics of the particle production, 
we treat the scalar QED under constant electric fields which is independent of r\ and neglect the 
back-reaction of pair creation to the background field. We drop the dependence of the electric 
fields on r and xj_ for simplicity. Under the constant electric field, we can derive the analytic 
solution for the distribution of created particles. Extension to the more general case, quark pair 
creation under a color-electric field and taking the back-reaction into account, is possible in a 
similar way as that in Ref. [13j . 

Pair production in electric fields confined to a bounded region has been studied by Martin 
and Vautherin [15] . They calculated the pair creation rate by using the Balian-Bloch expansion 
of Green's functions and concluded that if the boundaries move with the speed of light, finite-size 
correlations vanish and the pair creation rate is the same as Schwinger's original one. In their 
calculation, information from outside the light cone is not used. However, it should be used 
as an initial condition on the forward light cone because inside and outside the light cone are 
causally connected. Furthermore, calculating the pair creation rate is not sufficient to reveal 
the time evolution of the system [12]. To properly describe the system evolution, we must set 
an initial condition and investigate the time evolution of the quantum fields. A study on this 
line has been carried out by Cooper et al. [HI [IT]. Because the electric field is absent outside 
the forward light cone, there is no particle created. Therefore, the initial state on the light cone 
should be a vacuum. In Ref. [16] . the initial vacuum is defined so that a charge current, which 
is regularized by the adiabatic regularization, vanishes at an initial time. 

However, in quantum field theory, there is an ambiguity in identifying a vacuum. Requiring 
the absence of a particle or current at an initial time is not sufficient to determine one initial 
vacuum. Because we work in the flat Minkowski space, the vacuum must hold the Poincare 
invariance, even if we use the curvilinear coordinates. Hence, the vacuum we should employ for 
the particle mode defined in the r-rj coordinates must be the same as that associated with the 
usual quantization of free fields in the Cartesian coordinates. In other words, because particles 
captured by detectors are considered as the particle modes associated with the usual plane wave 
expansion of the field operator in the Cartesian coordinates, the state before the collision should 
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be the vacuum for that particle mode. 

Because the configuration of the electric field is not symmetric under translation along the 
beam direction (z axis) but symmetric under translation of the space-time rapidity rj, particles 
created from the electric field are eigenstates of the momentum conjugate to r/ (which will be 
denoted as A in this paper), instead of the usual longitudinal momentum p z . In Ref. [16], 
the relation between these two kinds of momentum variables is given based on a coordinate 
transformation in classical mechanics, and the rapidity independence of particle spectra is shown. 
In contrast, we will derive a relation between the two eigenstates of these momenta as an operator 
relation in a field-theoretical sense and show that the eigenstate of the momentum conjugate 
to T) is a quantum-mechanical superposition of the eigenstate of the longitudinal momentum 
p z . This prescription enables us to study multiparticle correlations. Because particles are 
created as a coherent superposition of several momentum modes, they are correlated in the 
momentum (rapidity) space. Two-particle correlations between particles produced in the early 
stage of heavy-ion collisions [18] attract interests in the context of the near-side ridge phenomena 
observed in heavy-ion collisions at the Relativistic Heavy Ion Collider (RHIC) [19} [20} \2l\ and 
recently in proton-proton collisions at the LHC [22]. Multiparticle correlations arising from the 
Schwinger mechanism in spatially homogeneous electric fields have been studied by Fukushima, 
Gelis and Lappi [23J. We will investigate how that correlation is modified if an electric field 
spans only inside the forward light cone. 

This paper is organized as follows. In Sec. El we review the field quantization in terms 
of the r-rj coordinates. Discussions in this section are to a large degree based on Ref. [14] . 
Because of the translational invariance to the r/ direction, particle modes are parametrized by 
the momentum conjugate to r]. In Sec. [3j the physical meaning of the momentum conjugate to 7] 
is investigated. After these provisions, in Sec. [H we will study particle production in the boost- 
invariant electric field which is confined in the forward light cone. The proper-time evolution 
of the momentum distribution functions is displayed there. Finally, two-particle correlations 
between particles created in the boost-invariant electric field are investigated in Sec. [5j 



2 Quantum fields in the t-tj coordinates 

In this section, we will review the canonical quantization in terms of the curvilinear t-ij coordi- 
nates, as a basis to study the pair creation in a boost-invariant electric field which spans only 
between two charged plates receding from each other at the speed of light. 
The t-t) coordinates (r, rj, Xj_) are related with the Cartesian coordinates as 



T = \Jt 2 - z 2 , r) = - in yj—;) ' Xj - = ( x ' y )- W 

In terms of the r-rj coordinates, the Klein-Gordon equation is written as 
/ d 2 d d 2 d 2 \ 

(j) is a charged scalar field with a mass m(^ 0) and a charge e. We quantize the field <j) on an 
equal r-line instead of an equal i-line, imposing the canonical commutation relation 

[0(r,r/,xj_),7r(r,r/',x|)] = i5(rj - r]')5 2 (x ± - Xj_), (3) 
where 7r(r, t],x±) = r^^(r, r], xj_) is canonical conjugate momentum. 
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Because of the translational invariance to the xj_ and rj directions, a solution of Eq. ([2]) can 
be expanded by Fourier modes as 



r,ri,x±)= / dXd p ± Pj _,a(t) 
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(4) 



This equation tells us that a particle mode defined on an equal r line can be parametrized by 
quantum numbers and A. To define this particle mode, let us decompose the field operator 
into a "positive frequency" part and a "negative frequency" part: 



dXd p± 



dXd p± 



Xp, ,x(r)a P± ,x + x 



«>,A( T ) b -Px,-A 



1 



where 



1 



. e «Px-x ±e «Ar7 



(5) 



(6) 



are c- number solutions of the field equation ([2]). The mode functions A (r, 77, xj_) are set to 
obey the orthonormal conditions 



J(A-A')5 2 ( P± -pi), 
. 



The inner product (</>i,</>2)r is defined as 



(7) 



>i, 02)r = i / dr]d 2 x ± T ( 4>\—4>2 ] , 



where 



a . ix _ A ,t 



'2 - ^r<Pi • <?2 



(8) 



(9) 



Owing to the orthonormal conditions ([7]), the commutation relations 



[ a p±,A,aL iA ,] 



5(A-A0^(p±-pi), 
[6px.A,bL A ,]=<5(A-A')5 2 (p±- P i), 



Pi, A' 
t a Px,A ; &p' ,A'] 



(10) 



[ a Px,A' b P',v] = ' 



are deduced from the canonical commutation relation ([3]), and the operators ap^ A and A 
acquire the role of the creation operator of a particle with quantum numbers pj_ and A. Fur- 
thermore, the expression of the charge operator in terms of a p±i A and b PXiA 



Q = e( 



s 



e / dXd p± 



a P± ,A a Px,A - bpx,A& Px ,A 



(11) 



indicates that A is the creation operator of a particle with charge +e and bj^ A is the creation 
operator of an antiparticle with charge — e. 
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In order that Eq. ([5]) is a solution of the Klein-Gordon equation ([2]), Xp ± \ ( T ) must satisfy 



(12) 



where m± is the transverse mass defined by m]_ = m 2 + pj_ . Note that the role of transverse 
degrees of freedom is only shifting the mass. Solutions of Eq. (|12|) can be expressed by the Bessel 
functions. For example, we may use the Bessel function of the first kind: 



X p± ,x(r) ~- 



(13) 



2sinh7r|A| ± ^ 

The normalization factor has been determined by the condition ([7]) or equivalently 



• I +t ___ + 

• f +t 1 - N 



1, 



(14) 



0. 



One can construct a Fock space by using the creation and annihilation operator associated with 
the solutions (|13p . However, the choice of solutions Xp ± \ is n °t unique. There remains the 
freedom of the Bogoliubov transformation: 

Xp ± ,x( T ) = «px,a Xp ± , A ( r ) + % ± ,\ X p± , A ( r ) 

X p± ,a( T ) = «Px,A Xp ±iA ( T ) + Xp ±>X ( T ), 

where the coefficients « p± ,a and /3 p± ,a satisfy the normalization condition 



(15) 



\ a p±,X\ 



Px,A| 



1. 



(16) 



For any set of such coefficients, Xp ± \ are a l so solutions of Eq. (|12p fulfilling the orthonormal 
condition (|14p . Therefore, there are infinite numbers of ways to decompose the field operator <j) 
into a "positive frequency" part and a "negative frequency" part, and each decomposition gives 
different sets of creation and annihilation operators: 



<Kt,??,X ± ) = / d\d 2 P± [xp ± ,A( r ) a P±,A + Xp ±i A( T )^- p± ,-A 



dXd 2 p± Xp ± , x ( T K±,x + X P± ,x( T )V- p± ,-x 



1 



1 



.giPi-XLg^Ar? 



.giPx-XXg^Ar? 



(17) 



The two kinds of creation and annihilation operators introduced above are related by the fol- 
lowing Bogoliubov transformation: 



dpx,A = a P x,Afi P x,A + P P j_,\b p±x 
t _ „,* ut 



b Px,A = °W Ax,A 



+ /3; ± ,Aa P x,A- 



(18) 
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Because creation and annihilation operators are mixed by the Bogoliubov transformation, the 
vacuum |0) defined by a PXi A|0) = b p± ,A|0) = and the vacuum |0) defined by 5 PX) a|0) = 
bp ± a|0) =0 are inequivalent: 

|0>^|6>. (19) 

That is, there are many inequivalent vacua in quantum field theory. This problem is not owing 
to the use of the special coordinates, the t-t] coordinates, but is inherent in quantum field theory. 
In the Cartesian coordinates, however, positive and negative frequency solutions are selected as 

<f£(t,x) = —= 1 e^rWp-x (20) 

according to the symmetry of the space and time. Now, oj p is a one-particle energy: oj p = 
\/p 2 + m 2 . In the t-T] coordinates, such a selection of positive and negative frequency solutions 
is not obvious because the metric depends on r and there is no translational symmetry on t. 

A proper set of positive and negative frequency solutions, in other words, a proper definition 
of particle in the r-r/ coordinates can be found in the following way [14J. The use of the r-rj 
coordinates is only a matter of description and should not change the physics. Therefore, even 
if one uses the r-ry coordinates, a physical vacuum must be the same as that defined in the 
Cartesian coordinates. That is, a "positive frequency" solution in the r-rj coordinates must be a 
superposition of only the positive frequency solutions defined in the Cartesian coordinates; the 
negative frequency cannot be mixed. We can construct the solutions satisfying such a condition 
by using the Hankel functions as follows: 

e ip±-x ± +iXri 

<PLx(t,V,x±) = X± )A (t)^==-, (21) 

where 



X+, A (t) = ^ef (22) 
X- .» = \xt ,(r)l * = -^§e^hm ± r). (23) 



2? 



Using the integral representation of the Hankel function 

j p niv/2 poo 

HP(z) = / dt e -^^ s ht+ut y Re ^| < l lmz < j j ^ 



7T 

L± 



one can rewrite (j) A (r, rj,x.±) as follows: 



^■sjC^**' <26) 



where y p denotes rapidity corresponding to momentum p: y p = tanh.^ 1 (p z /ui p ). The condition 
Im z < in (124p is satisfied if we suppose that the mass has a small imaginary part asm - ie, 
which is a usual prescription in quantum field theory. An important point is that the positive 
frequency and the negative frequency solutions are not mixed with each other. In the forward 
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light cone, one can decompose the field operator in two ways by using the mode solutions in the 
Cartesian coordinates p (i, x) [Eq. (j20l) ] and in the t-ij coordinates <j> a (t, r/, xj_) [Eq. pT) ]: 



0+(i,x)ct p + p (i,x)6L p 



(r,r?,x_ L )a px 



>,77,x ± )bi p±i _ A 



(27) 
(28) 



Then, the relations between the creation or annihilation operator associated with </>p(i,x) and 



that with 



'p±,a 



V, r],x.±) are obtained as 



1 



hi 



Their inverse relations are 



1 



dX e iXy "a- 



p ± ,A i 



1 



dp, ,A 



P±,A 



(29) 



P±,A 



2vr 
1 



-*Ayp 



oo V w p 

OO 



27T J-oo v^P 



(30) 



From these equations, we see that the particle mode respecting the boost invariance is a super- 
position of various momentum modes. Because the particle modes and the antiparticle modes 
are not mixed with each other, a vacuum defined in the Cartesian coordinates a p |0) = 6 P |0) = 
is also a vacuum for the particle mode defined in the r-rj coordinates: 



a Px ,A|0) = b PXiA |0) = 0. 



(31) 

Using these particle modes defined by Eq. ([29]) . we can properly calculate field-theoretical quan- 
tities in terms of the T-rj coordinates. 



3 The meaning of the momentum conjugate to the space-time 
rapidity 

In the previous section, we have formulated the field quantization in terms of the r-rj coordinates 
and derived the relation between the particle modes respecting the boost invariance and those 
associated with the usual plane wave [Eq. (|29p ]. In that process, the quantum number A has 
been introduced as a conjugate variable to r/. However, the physical meaning of A is not obvious 
from the discussions in the previous section. We will investigate it in this section. 

3.1 Classical mechanics 

First, we study the meaning of A in the framework of the classical mechanics in the r-rj coor- 
dinates. For simplicity, we analyze a (l+l)-dimensional system. The action of a free particle 
is 

S = -m / dty/l - z 2 = -m dT\/l-T 2 f] 2 , (32) 
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where 

dz , , 

"* (33) 

* = £ < M > 

Hence, the Lagrangian in terms of the t-t] coordinates is 

L T = — my/l — T 2 f] 2 . (35) 

Motion of the particle is expressed as z = z(t) or rj = f/(r), both of which give the same motion 
if the particle is in the forward light cone. The momentum conjugate to 77 is 

X = ^= ^ ■ (36) 



Using the relation 

one can rewrite Eq. (|36p as 



V = —. 7, 37) 

r t — zz 

m(zt — z) 

A = — = p 2 t - w p 2:, (38) 
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where p 2 is the momentum conjugate to z and uj p is the energy of the particle: 

mz , , 

^ " (39) 

(40) 



To make the meaning of A clearer, let us boost this system by the velocity v z = z(t)/t in 
the z direction, where z(t) is the position of the particle at time t. By this boost, t and z are 
transformed as 

t^t' = t (41) 

z -»• z' = 0, (42) 

where the prime denotes a boosted quantity. Because A is a boost-invariant quantity, 

A = p z t — ujpZ = p z t' — uj'pZ 1 = p' z T . (43) 

This equation clearly tells us the meaning of A: A/r is the longitudinal momentum of the 
particle in the frame which moves with the longitudinal velocity v z = z(t)/t. Note that the 
velocity v z = z(t)/t is time-dependent unless the motion of the particle is a straight line from 
the origin z(t) oc t. That is why p' z = A/t depends on r even if the particle is free and A is 
constant. 

Under the gauge field A T = 0, A v = \Et 2 giving a constant electric field, the action of a 
charged particle is 



S = —m 

which leads the equation of motion 



J drs/l - r 2 i] 2 - e J dr A v rj , (44) 

-(A-i^.O. (45) 

If A = at r = 0, this equation implies that p' z = A/r = ^eEr. It seems that the acceleration 
the particle feels is half compared with the description in terms of the Cartesian coordinates: 
p z = eEt. This is because v z = z(t)/t varies under the electric field and thereby the frame where 
the longitudinal momentum p' z = A/r is defined is also changing in time. 
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3.2 Quantum field theory 

In this subsection, we investigate the physical meaning of A in the frame work of quantum field 
theory. The key of this analysis is Eq. (|29p . which relate the particle modes in terms of the r-rj 
coordinates and those in the Cartesian coordinates. 

Let us define a one-particle state with quantum number A and qj_ as 
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|qi,A) = y^c4 iA |0>, (46) 

where the normalization factor is introduced so that (qj_, A|qj_, A) = 1, and is a space volume 
on the r-constant hypersurface, introduced as 

V n = L V L 2 = J rdi] J d 2 x± 

= t Jdr] e^-^ jd 2 x ± ^(px-px)-^ ( 47 ) 

= r(2vr) 3 <5(A-A) ( 5 2 (p ± -p ± ). 

By calculating expectation values of several quantities with this one-particle state, we shall 
examine the meaning of A. 

First, we study the expectations of particle number operators. The expectation of the number 
operator in the a p± ,\ basis is 

dN I \\ t I u 

^V=^' A l a Px,A'W|ciL,A) 

= ^[*(A-A> 2 (p,-qx)] 2 <«> 
= 5(A-A , )<5 2 (p±-qi). 

This is a trivial result since the expectation is taken by the eigenstate of pj_ and A, and does 
not tell us anything about the physical meaning of the quantum number A. 
Next, we study the expectation of the number operator in the a q -basis: 

dN I \ I t I x\ 
-3- = (qj_, A|a p a p |qj_, A) 



d?p 

= 2^"/ dX '^ X ' Vp J dAV^q^Ala^ap^Hq^A) 

1 r(2vr) 2 rr2 , .2 
UJp v v 
1 T 



(49) 



-— 5 2 (p ± -c LL ). 

Converting the longitudinal momentum p z to the rapidity y p , we can get the rapidity distribution 

dN _ dN 

dy p d 2 p ± p d 3 p , 5Q s 

= ^ 2 (p ± -qJ. 

v 

This is independent of rapidity y p . The one-particle state |qj_,A) contains all rapidity modes 
with equal weight. This result is consistent with the boost invariance of the state |qj_, A) in the 
longitudinal direction. 
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Furthermore, because several momentum modes are condensed in the state |qj_,A), 



(q±, A|a+a p /|qi, A) 



V, 



(51) 



is nonzero even if p z ^ p' z . This implies that nontrivial correlations in the longitudinal momen- 
tum space may arise from the particle production in the boost-invariant electric field. We will 
investigate it in Sec. [5j 

Because the rapidity distribution (|50p is also independent of A, the meaning of A is not 
obtained from it. To further investigate the meaning of A, we study the energy- momentum 
tensor. In the Cartesian coordinates, the energy-momentum tensor for free complex scalar fields 
is 



-g 



(52) 



(/i, v = 0, 1, 2, 3), where g ta> = diag(l, —1, —1, —1) is the metric of the flat Cartesian coordinates 
system and \J—g = y^detg^l = 1. Also in the T-rj coordinates, it has the same form: 



-0 



L 



+ 



m 



(53) 



(a,/3 = r, 1 , 2, 77) , where Q a ^ = diag(l, — 1, — 1, — 1/t 2 ) is the metric of the T-rj coordinates 
system, and y/^ 



9 = y/\detg a p\ = r. 

The relations between the differentials in the Cartesian coordinates and those in the T-rj 
coordinates are 



(54) 



/ cosh rj — sinh ?A 

10 

1 

\— sinh 77 cosh 7/ / 

Note that the matrix in the right hand side accords with the matrix representing the Lorentz 
boost along the z axis with the velocity v z = z/t = tanhr/ : 



d 1 
d 2 



/d°\ 
d 1 
d 2 

w 



Therefore, 





/ cosh rj 





— sinh rj\ 











1 














1 









V- 


- sinh rj 





coshr/ J 










rprl 




1 






rpll 


T 12 tT 1t > 




3 


rp2r 


rp21 


T 22 T rp2r] 








rT r,l 


tT^ 2 t 2 


jyn J 



(55) 



(56) 



equals the energy-momentum tensor transformed by the Lorentz boost 

fi* = A^A u p T ap . 



(57) 



For later reference, let us define the rj- frame as the frame which is boosted by (J55J). Then, 
T^(t, r/) is the energy- momentum tensor observed in the rj- frame. 

The expectation value of T 03 (r, rj) = T TV (T,rj), which represents the momentum density in 
the z direction observed in the 77-frame, with the one-particle state |pj_,A) is 



(p±,X\T T \T,rj)\ P± ,X) 



A 1 



(58) 
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This means that the momentum the state |pj_, A) contains is A/t. This result is consistent with 
the result of the classical mechanics in the previous subsection. In the present case, however, 
7] is a parameter indicating a space point, while in the classical mechanics 7/ is a mechanical 
variable representing a point where a particle exists. Therefore, Eq. (I58p means that momentum 
A/t is distributed at any space point of r\. 



4 Pair creation in the forward light cone 

In this section, we study the pair creation in a constant electric field which exists only inside 
the forward light cone. 

The metric of the t-tj coordinates has a singularity in r = 0. To avoid this singularity, we 
suppose that the electric field is zero when < r < To and is switched on at r = to > 0. Such 
an electric field is given by the gauge 

\i Er ( t > t o) 

and A T = 0. The pure gauge A v = \Etq at t < to is introduced so that the gauge field is 
continuous. It will be confirmed later that if To is sufficiently small, results of the calculations 
are insensitive to the values of To- 

The Klein-Gordon equation under this gauge field is 

[t 2 8 2 + t8 t - D 2 - T 2 dl + m 2 r 2 ] <f>(r,rj,x ± ) = 0, (60) 

where D v = 8^ + ieA^. A solution of this equation can be expanded in the same way as the free 
field case: 

0(t, rj, xl) = J d 2 p ± d\ [0+ in A (T, 77, x ± )< ±iA + ^ a (t, v, Xl)^J ±t _ x ] . (61) 

(f > p ™x( T j V? X J-) i s a positive or negative frequency mode function satisfying Eq. (|60|) . The su- 
perscripts 'in' specify the initial condition for the field: At t < to, there is no electric field so 
that the mode functions are free ones. The free solutions at t < to are given by the gauge 
transformation (A v = -> Et^/2) of Eq. (j2T|) : 

V V ; (62) 

The mode functions 4>%™\ at t > To are constructed so that pj "\ and their derivative with 
respect to r are continuous, respectively, at r = tq. Their explicit forms can be found in the 
appendix [A] By expanding the field operator with these mode functions <ftp*\ as Eq. (|6"Tj) . we 
can obtain the annihilation operators op^ A and b 1 ^ A and the associated vacuum state |0, in), 
which is defined as 

a^ iA |0,in) = b^ iA |0,in) = 0. (63) 

We choose this vacuum as the state of this system. This choice is a manifestation of the initial 
condition that there is no particle before the electric field is turned on. An important point 
is that this vacuum is also the vacuum for particles defined in the Cartesian coordinates as 
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discussed in Sec. [2j Therefore, although the r-77 coordinates can cover only inside the forward 
light cone, the absence of particles outside the forward light cone is guaranteed. 

Because we now use the Heisenberg representation, the state |0, in) does not change during 
time evolution. What evolves in time is a definition of particles or, in other words, creation and 
annihilation operators of particles. Such a time-dependent definition of particles is introduced 
by decomposing the field operator into positive and negative frequency instantaneously [12]: 



d p±d\ 



^S ) (T,i7,x J .)a BL ,A(n) + ^S^r,f / ,x J .)b^ if _ A (r l ; 



(64) 



where 



? P±,A 



(r, 77, xj_) is a positive or negative frequency solution of the equation of motion ([SO 



under the pure gauge A v = A v (t = ti): 
< ± ( ,A } (W ± ) 



1 



2i 7(2^)3 



e *P± -xj. g i[A-eA^ (ri )]r? 



(ti), 



1 



gip±-xx e «[A-eA^(ri)]r;^ 



(65) 



The operators CI p±) a(t) and b p±i \(r) give a time-dependent particle definition. The instantaneous 
mode functions (r, ry,xj_) are related with the in mode functions 4>p^\(T,rj,x.±) as follows: 



(t^xjJ 



-(n) 



Pi. 



Px,A 



; ap J _,A+eA^(ri)-eA^(r )(Tl) ^px,A+eA^(Ti)-eA^(To) t T > 
+ / 5 Px,A+eA ?) (Ti)-eA ?) (ro)( r l) ^ Px ,A+eA^(-n)-eA,(r ) ( T ' "7> Xj -) 



Px 



(66) 



pi.A'^Px^ 



9 p|,A' ' V,A 



/?px , A+eA„ (ri ) -eA„ (to ) ( r l ) ^px^A+eA, (n )-eA„ (r ) ( T ' ^ ' Xj - ) 

v T l) \i»4„fT^-R/l.,f'-m , >( T > ? ?> X x, 



+ a 



Px,A+eA,,(ri)-eA^(ro) 



where the Bogoliubov coefficients a p±i A( r ) and (^p ± .x{ T ) are introduced bvRI 



(67) 



t+(n) - + ir 

Vi.A" V,a; t 

/-in 

V ,A" V,A 



api,A'(n)5 2 (p± - pi)5(A - A' + eA„(ri) - e^(r )), 
/3 p ;,A'(nM 2 ( P x - pi)5(A - A' + e^(n) - eA v (r )). 



(68) 



For the explicit forms of the Bogoliubov coefficients, see Eqs. (|A18|) and (|A19p . These coefficients 
satisfy the condition |ct Px ,A | 2 -|/?px,a| 2 = I- By the definition of ^(r), /3 Px A (r) = for r < r . 

Inserting Eqs. (|66p and (|67p into Eq. (|6ip and comparing it with Eq. (|64p . we get the rela- 
tion between the particle definition associated with in-solutions and the instantaneous particle 
definition: 



0px,a(t) - apx,A(T)ap n L ,A-eA^(r)+eA^(ro) + Pp±A T ) b -J ± -X+eA v (T)-eA v 



(r )' 



-P-L.-AV / 



o 



Px,A 



V1b in f 

^ > -Px ,-A+eA^(r)- e A^(ro) 



+ /3 P x,AWa: 



(69) 



Px,A-eA,,(r)+eA^(T )' 



* 1 Actually, we should write otp ± ,\{j, Ti) instead of q Pj _.a(ti) in Eq. (|68|) . However, we omit one time-argument 
for simplicity, because the instantaneous decomposition of the field [Eq. (|64|l ] is physically relevant only at r = ti. 
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Figure 1: (color online). Longitudinal momentum distributions with fixed transverse momentum 
(V^Et = 0.1) 




Figure 2: (color online). Transverse momentum distribution with fixed longitudinal momentum 
X/VeEr = 1 (m 2 /2eE = 0.1 and y/eEr = 0.1) 



Because the creation and the annihilation operators are mixed by the Bogoliubov transfor- 
mation (I69p . the vacuum expectation of the number operator can be nonzero: 

^ = <°.H"L,aMwMI0>> 

This means that particles are created from the electric field. Because of the charge and the 
momentum conservation, antiparticles have always opposite momentum to particles: 

(0,in|b^ A (r)b p±)A (r)|0,in) = \(3_ p±i _ x (r)\ 2 (71) 

In Figs. 11121 the proper-time evolution of the distribution function 



is plotted. Figure Q] shows the longitudinal momentum distributions with fixed transverse mo- 
mentum pj_ = for m 2 /2eE = 0.1 and 0.01. The longitudinal momentum distributions are 
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plotted as a function of A/r instead of A because A/r denotes a momentum in the rj- frame as 
explained in Sec. El Figure [2] exhibits the transverse momentum distribution with fixed lon- 
gitudinal momentum X/y/eEr = 1 for m 2 /2eE = 0.1. Hereafter, all figures are shown in a 
dimensionless unit scaled by \J eE. 

These momentum distributions are very similar to those in the uniform electric fields (see 
Ref. p]): 

• Particles are created with approximately zero longitudinal momentum. 



• In the transverse direction, the distributions are nearly Gaussian exp 




• After being created, they are accelerated to the longitudinal direction by the electric field 
according to the classical equation of motion. 

The similarity of the transverse momentum distributions is quite reasonable, because the fields 
are free in the transverse direction also in the present case. In contrast, the similarity of the 
longitudinal momentum distributions is only a superficial one. The physical contents are quite 
different between the uniform electric field case and the present case, because the physical 
meanings of the longitudinal momentum are distinct. In the uniform case, the longitudinal 
momentum is defined in the center of mass frame. On the other hand, in the present case, the 
longitudinal momentum A/r denotes the momentum observed in the 77- frame, which is the frame 
boosted by (|55p from the center of mass frame. Therefore, Fig. [1] does not mean particles are 
created with zero longitudinal momentum in the center of mass frame, but it means particles 
are created with the scaling velocity distributions. That is, a particle created at the point 
(i, z) has the velocity z/t from the first instance when it is created. This result supports the 
assumption in Refs. OHO] that a source term in a kinetic equation contains the factor 5{y — rj). 
In our field-theoretical treatment, however, the velocity is not exactly v z = z/t because of 
quantum fluctuation. After being created, particles undergo the acceleration by the field and 
their velocity distribution deviates from the scaling one, of which processes are expressed by 
A = \eEr 2 + const, in a classical level. 

The momentum distribution with respect to momentum defined in the Cartesian coordinates, 
in other words, momentum in the center of mass frame, is calculated as follows: 

^ = (0, in| aj,a p 1 0, in) 

= 2hrJ dX I dX '^ t{ ^ X ' )y H0M^(r)a p±i y(T)\0M) (73) 

= 2^ l dX /p - A(T) (27F' 
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Figure 4: (color online). Dependence on the time of the switch-on tq (a = 0.1) 



in which derivation Eq. (j29j) was used. Thus, the rapidity distribution 

d^ P± dy p = = Wfl /P± ' A(T) {7) 

is independent of rapidity. Because we have assumed the perfect boost invariance, the central 
plateau in the rapidity distribution extends to infinity. In reality, the speed of nuclei after 
a collision is not exactly the speed of light, so that the configuration of the electric field is 
not perfectly boost-invariant, and there should be some cutoff in the plateau. Note that this 
independence of rapidity is a direct consequence of Eq. (j29|) and is not affected by the explicit 
form of / p±! a(t). Therefore, even if the electric field has proper-time dependence or the back- 
reaction is taken into account, the rapidity distribution (I74D is always independent of rapidity 
as long as the boost invariance of the electric field is assumed. 

In Fig. [3J the time evolution of the rapidity distributions ([711) divided by r is plotted. They 
show linear increase at later time. That is, the particle number density increases quadratically in 
proper time. This is because (i) pair creation happens constantly some time after the switch-on 
of the field and (ii) the space volume the electric field spans increases linearly in t. 

The dependence on the time of the switch-on to is shown in Fig. [U Figure H^a) represents 
the longitudinal momentum distribution at fixed time y/eEr = 10, and Fig. H{b) does the time 
evolution of the particle number density. Three lines corresponding to tq = 0.01,0.1 and 0.5 
are almost overlapped in both figures. These demonstrate insensitivity to the values of To for 
sufficiently small tq (V eEro < 0.1). 



5 Two-particle correlations 

Because particles are created as a coherent superposition of several momentum modes in the 
electric field which exists only inside the forward light cone, they are correlated in the momentum 
space. In this section, we study the two-particle correlations in the momentum space between 
particles created from the boost-invariant field. 
The two-particle spectrum is defined as 

(0, inja^aJ-apaqlO, in). (75) 



d 3 pd 3 q 



In the case of the pair creation in a spatially uniform electric field [23], the two-particle spectrum 
is 

dN 2 dNdN 2 V x3 



d?pd?q d 3 pd 3 q n (2vr) 



+ {/p(*)}^(P- ( l). ( 76 ) 
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where f p (t) is the momentum distribution function of particles created from the field. The 
first term in the right-hand side is an uncorrelated part, which is a product of one-particle 
distribution ^ = (0, in|apa p |0, in) = f p (t) , and the second term is a correlated part. 
Particles are correlated only if they have identical momenta. This correlation is due to the 
Bose-Einstein statistics. 

Let us see how this two-particle spectrum is modified if the electric field localizes in the 
forward light cone. Using Eqs. (|29p and (|69p . one can derive 



_dN^_ 

(PpdPq 



1 



dN dN 



d?p d 3 q ujpUJg 



dX 
2~77 



e -^-^)/ pXjA ( r ) 



L 2 



r<5 2 (p± - qj_) 



instead of Eq. (|76|) . Converting the longitudinal momenta to rapidities gives 



(77) 



dN? 



dN 



dN 



d 2 p x dy„d 2 qi_ dy q d 2 p 1 _ dy„ d 2 q±_ dy c 



+ 



d\ 

277 ( 



-i\(y p -y q ) 



/px,a(t) 



L 2 



r5 2 (p±-qx). (78) 



For later reference, we introduce the correlation function and its longitudinal part as 



C(yp,p±;y<?,qL) 



dN 2 

d 2 p±dy p d 2 q ± dy q 



dN dN 
d 2 P±dy p d 2 q ± dy q 



dN dN 
d 2 p ± dy p d 2 q L dy q 



\f£ e -*Hy P -y q ) fp±AT) \ 2 6 2 {p± 



(79) 



C L (Ay = y p - y q ,p±) 



Ufe-^/pxAWl 



(80) 



{/f|/px,A(r)} Z 

From Eq. (|77p we can see that the correlation in the transverse direction is the same as that in 
qj_). This correspondence is reasonable because in the transverse direction 
This correlation is short-range in momentum space. What is 



Eq. (EH: j^ 7 6 2 ( P± 
the field is free in both cases 



remarkable is the correlation in the longitudinal direction. The correlated part no longer contains 
the delta function 5{p z — q z ) and is given by the Fourier transform of the momentum distribution 
. Jd± e -iX(y p -y q )f^ x ( T y Therefore, particles are correlated even if their longitudinal momenta 



are different, while in Eq. (176)) particles are correlated only if they have an identical momentum. 
Because the rapidity correlation is given by the Fourier transform of the momentum distribution, 
its width is approximately the inverse of the width of momentum distribution f p±t x(j) in A space. 
In particular, if the distribution is proportional to S(X), in other words, if the particles' velocity 
distribution is the scaling one, C^(Ay, pj_) = 1 for any Ay and the rapidity-correlation range is 
infinite. 

Figure [5] exhibits the time evolution of (a) the longitudinal correlation Cl under the constant 
electric field and (b) the corresponding momentum distributions f p±! \, whose Fourier transfor- 
mations give the correlation functions through Eq. (180|) , Because of the pair creation and 
subsequent acceleration, the width of the longitudinal momentum distribution grows in time. 
As a result, the longitudinal correlation fades away as time goes on. This can be understood 
as follows: Because particles which are created at the same time are strongly correlated, the 
correlation is long-range at first. As particle production continues to happen, many particles 
which are created at different times begin to coexist, so that the correlations are attenuated at 
later times. 

2 

In Fig. El the a = dependence of (a) the longitudinal correlations at fixed time V eEr = 1 
and (b) the corresponding momentum distributions are shown. The two peak structures seen 
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(a) the longitudinal rapidity correlations Ci(Ay) (b) the longitudinal distributions / p±) a(t) 

Figure 5: (color online). Time dependence of the correlations and the corresponding momentum 
distributions (m\/2eE = 0.1 and \feE~TQ = 0.1) 




(a) the longitudinal rapidity correlations Ci(Ay) (b) the momentum distributions f PXl x( T ) 

2 

Figure 6: (color online), a = t^jj dependence of the correlations and the corresponding momen- 
tum distributions (V eEr = 1 and y/eEro = 0.1) 



in the momentum distributions at early time are owing to some quantum effects and have been 
obtained in earlier works treating the pair creation of scalar particles |12t 123] . The rapidity- 
correlation length decreases with increasing a. This is because the height of the momentum 
distributions depends on a exponentially (~ e~ 2na ) [12] and their width scaled by V 'eE is nearly 
independent of a, so that the relative width of the momentum distributions increases with 
increasing a. This a dependence means that the heavier the transverse mass or the weaker the 
electric field, the shorter the rapidity correlation. 

Before closing this section, let us note that although we have dealt with the constant electric 
field which is independent of r, the formula (|77p is valid for electric fields having arbitrary r 
dependence. Thus, we can use Eq. (I77j) even if the back-reaction of the particle production 
is taken into account, within the approximation that interaction between produced particles 
through quantum gauge fields is neglected. 



6 Summary and discussion 

In this paper, we have studied the pair creation of scalar particles in electric fields which span 
only inside the forward light cone as a model for plasma formation in heavy- ion collisions. We 
have assumed the electric field is symmetric under the Lorentz boost along the longitudinal beam 
direction and investigated how particle production happens maintaining the boost invariance of 
the field. 

To give a description which holds the boost invariance, we have developed the field quanti- 
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zation in terms of the curvilinear t-t] coordinates. This formalism may be useful for not only 
the study of pair creation but also for other studies that quantum fields in a boost-invariantly 
expanding system are involved. Although we have treated scalar fields for simplicity, general- 
ization to Dirac fields or vector fields is straightforward. A detailed formulation of Dirac fields 
in the T-r\ coordinates has been presented in Ref. [25] . 

We have shown that under the boost-invariant field, particles are created as a coherent 
superposition of several eigenmodes of the longitudinal momentum. As a result, the rapidity 
distribution of the particles is independent of rapidity, and the particles have the scaling velocity 
distribution v z = z/t from the first instance they are created. A significant point is that this 
flow of particles consists of quantum-mechanical states. Each particle spreads in the forward 
light cone and forms the scaling velocity distribution as a quantum state. It is not the case that 
a cluster of classical particles forms this flow like in the Bjorken flow in classical hydrodynamics. 
This fact brings nontrivial rapidity correlations among the particles. 

We have calculated the two-particle correlations between particles created from the boost- 
invariant field, and found that the correlation is short-range with respect to the transverse 
momentum, which originates in the Bose-Einstein correlation, and is long-range with respect 
to the longitudinal rapidity. These features may remind us of the near-side ridge phenomena 
observed in nucleus-nucleus collisions at the RHIC [19j [20], [21] and recently in proton-proton 
collisions at the LHC [22J. However, to make quantitative predictions, further improvements of 
this study, such as introducing a proper-time dependence of the electric field due to expansion 
of the system and taking the back reaction of the pair creation into account, are necessary. The 
back reaction can be introduced by the Maxwell equation 



where j v denotes the charge current operator of the produced particles. In our calculation, the 
longitudinal correlation fades away as proper time goes on, because under the constant field the 
particle production continues to happen and many particles which are created at different times 
coexist at later times. Hence, we can conclude that the life time of the electric field should be 
short to obtain a long-range rapidity correlation. Furthermore, considering a tubelike structure 
of the electric field and effects of radial flow may be important for the transverse correlation |18j . 
It is also interesting to study whether the correlation can survive under a thermalization process, 
which is not included in our treatment, or under the subsequent hydrodynamic evolution. 
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Appendix A Solutions of the Klein— Gordon equation under 

the boost-invariant field 

In this appendix, we show a set of solutions of the Klein-Gordon equation [Eq. (|60p] 




(81) 



[r 2 d 2 T + rd T -D%- r 2 dl + roV] 0± iA (r, r,, x ± ) = 



(Al) 
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under the gauge A v = \Et 2 with the initial condition such that they are continuous with 



respectively, at r = To, where 



(0), 



2/ 



-(0)/ \ [ +(0)/ N 



(A2) 

(A3) 
(A4) 



Under the gauge potential = ^Et 2 , a set of normalized solutions of Eq. (|A1[) can be 



constructed by the confluent hypergeometric functions U (a, 6, z) [21] as follows: 



(t,7/,Xl) = Xp ± , A ( T ) 



(A5) 



where 



xp~ x , A ( r ) 



— A 

_L e -f («+A) 5 e -!^V(l/2 + ia + i\,l + iX, ieEr 2 /2) (A6) 



X p± ,a( t )= Xp ±jA W 



(A7) 



These solutions ^ A are not continuous with Eq. (|A2|) at r = To- To find solutions which are 
lily t( 

are defined by Eq. ([65]) : 



connected smoothly to Eq. (|A2p at r = tq, we decompose the solutions 0^ A by (j>p^x, which 



0+, a (t, 17, x ± ) = y d 2 ^A' [« i ( J ) ,^ x , A J r ^ T A V(r, ry, x ± ) - 
= V.A+eA^n) (r, n) 0p ± ( iA | eyli)(Tl) (r, r/, x ± ) 



7",77,X_Lj 



(A8) 



y Rj_,A 



(T,r/,xj_) 



pi.A-^^^^^X^)- 
V,A+e^(n)( T > T ^p ± ( , A +eA,(n)( T ' ^ X ^ 



,.-(ri) 



(n), 



T,J? ) XJ. / 



+ ^ P± ,A+e^(ri) 



(A9) 



where the Bogoliubov coefficients A p±j \(t,ti) and B pX; \(r, t%) are defined, respectively, by 



A p± , x (r,r 1 ) = l r{ X ^%) 



* d _ + 

RL,A-eA7(n) 

* ^ 



(r) 



(A10) 
(All) 



Substituting Eqs. (|A3p . (|A4j) and (jA6J) into these equations, one can derive the explicit forms of 
the Bogoliubov coefficients. Albeit Eqs. (|A8|) and (|A9P are mathematically valid for any positive 
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r and ti, the expansion by the mode functions 0p A (r, 77, xj_) is physically meaningful only at 
r = n. Therefore, it is sufficient to know the Bogoliubov coefficients having the same time 
argument r = t\\ A p±j \(t) = A p±t \(T, r) and -B PXj a(t) = B p±t \(T,T). After completing the 
differentiation in Eqs. (jAlOj) and (jAlip . we can set r = t\ and obtain the following equations: 



A p±>x (r) = ^e-^~^) (leEr^^K-^ 

x ^(l + 2ia)H W iX (m ± T)U(^ + ia + i(X-^eET 2 ),i(X-^eEr 2 ),^eET 2 ) (A12) 
-mrrH^im^U^ + ia + i(\ - \eEr 2 ), 1 + i(A - -e^r 2 ), ^ e ^ 2 )} 
^, A (r) = ^ e -f(«-|^ 2 ) (l e ^) t(A " |e " 2) 

x I (1 + 2ia)H$ (m ± T)U(^ + ia + i(X - ^eEr 2 ),i(X - ^eEr 2 ), ^Et 2 ) 
+m ± TH[% x (m ± T)U(^ + ia + i(X - \eEr 2 ), 1 + i(X - \eEr 2 ), l -eEr 2 ) 



(A13) 

Taking the relations (|A8j) and (|A9[) . we can construct the mode functions A which are 
continuous with Eq. (|A2p (which is equivalent to 4>^ T \) at r = to as follows: 

<Pp ± ,x(T,V,*±) = A * Px ,\(t )4>+ ± :A _ eA(To) (r,77,x ± ) - S; ±)A (T O )0p i)A _ eA(To) (r,77,x ± ) (A14) 

^P ± ,a( T '^ X ^) = A P±A T 0) 4> p± ,X-eA(r )( T ^^^ ~ B P± -A ( T o) 0+ , A _ e A(r ) ( T ' X ^ ) " ( A15 ) 

Further using Eqs. (fAHJ) and (|A9|) . we can expand the mode functions A by the instantaneous 
free solutions at r = t\: 

4> p± ,\( T ,V, x ±) = {- 4 p ± ,A( T 0)-4p ± ,A+eA^(r 1 )-e J 4 I) (r )(T"l) ~ -B P± ,A ( r o)-Bp ± ,A+e J 4 I) (ri )-eA t) (r ) ( T l ) } 
X ^L^+eA v {T 1 )-eA(T )^ T ^ r l^) 
+ { A P ± ,\( T o) B P± ,\+eA v (T 1 )-eA v (T )( T l) ~ B p±A T o) A P ±,X+eA v (T 1 )~eA v (r )( T l)} 

X ^piA+eA^rO-eACro)^' 7 ?'^) 
— a pL,A+eA^(ri)— eA,, (to ) ( r l ) <P p± ,\ +e A v ( n ) -eA(r ) ( T ' ^ Xj - ) 

+ / 3 P i ,A+eA t) (r 1 )- e A,(r )( T l)'/' p± ( i A+ e A I) (r 1 )- e A(r )( T ' »?> X ^) ( A16 ) 
^ p± ,a( T ' r ?' X J-) = a P x,A+eA 7 ,(ri)-eA ?) (ro)( T l)^p x ( ,A+eA I ,(ri)-eA(To)( T ' 7 ?' X ^) 

+ ^P±^+eA v (T 1 )~eA v (To)(n)4 l p ^x+eA v (T 1 )-eA{T )( T ^ ^ Xj - ), (A17) 
where the Bogoliubov coefficients 

a p± ,A(r) = A* X)A _ ej4?;(r)+ej47)(7 . o) (ro)Ap Xi A(T) - ^p^A-eA^+eA^ro^O^PL^CO (A18) 
/3 Pj _,a(t) = A p±iA _ ej 4 i) ( T ) +ej 4 i) ( TO )(To)-Bp ±iA (r) - B pj _ i X-eA ri (T)+eA rl (To){' r o) A P x,\( T ) ( A19 ) 
have been introduced. 
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